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§1. Introduction 

This paper contains a step toward a complete solution (yet to be achieved) of the 
celebrated Boltzmann-Sinai ergodic hypothesis. In a loose form, as attributed to L. 
Boltzmann back in the 1880's, it asserts that gases of hard balls are ergodic. In a 
precise form, which is due to Ya. Sinai in 1963 [Sin(1963)], it states that the gas of 
N > 2 identical hard balls (of small radius) on a torus T", v > 2, (a z/-dimensional 
box with periodic boundary conditions) is ergodic, provided that certain necessary 
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reductions have been made. The latter means that one fixes the total energy, sets 
the total momentum to zero, and restricts the center of mass to a certain discrete 
lattice within the torus. The assumption of a small radius is necessary to have the 
configuration space connected. 

Sinai himself pioneered rigorous mathematical studies of hard ball gases by prov- 
ing the hyperbolicity and ergodicity for the case N = 2 and v = 2 in his seminal 
paper [Sin(1970)], where he laid down the foundation of the modern theory of 
chaotic billiards. Then Chernov and Sinai extended this result to (N = 2, v > 2), 
as well as proved a general theorem on "local" ergodicity applicable to systems of 
N > 2 balls [S-Ch(1987)j; the latter became instrumental in the subsequent stud- 
ies. The case N > 2 is substantially more difficult than that of N = 2 because, 
while the system of two balls reduces to a billiard with strictly convex (spherical) 
boundary, which guarantees strong hyperbolicity, the gases of iV > 2 balls reduce to 
billiards with convex, but not strictly convex, boundary (the latter is a finite union 
of cylinders) - those are characterized by very weak hyperbolicity. 

Further development has been mostly due to A. Kramli, D. Szasz, and the present 
author. We proved hyperbolicity and ergodicity for A = 3 balls in any dimension 
[K-S-Sz(1991)] by employing the "local" ergodic theorem of Chernov and Sinai, 
and carefully analyzing all possible degeneracies in the dynamics to obtain "global" 
ergodicity. We extended our results to A = 4 balls in dimension v > 3 next year, 
and then I proved the ergodicity whenever N < u (this covers systems with an 
arbitrary number of balls, but only in spaces of high enough dimension, which is a 
restrictive condition). At this point, the existing methods could no longer handle 
any new cases, because the analysis of the degeneracies became overly complicated. 
It was clear that further progress should involve novel ideas. 

A breakthrough was achieved by Szasz and myself, when we employed the meth- 
ods of algebraic geometry [S-Sz(1999)]. We assumed that the balls had arbitrary 
masses mi, . . .,mjv (but the same radius r). Now by taking the limit — > 0, 
we were able to reduce the dynamics of A balls to the motion of A — 1 balls, 
thus utilizing a natural induction on A. Then algebro-geometric methods allowed 
us to effectively analyze all possible degeneracies, but only for typical (generic) 
vectors of "external" parameters (mi, . . . ,mjv, r); the latter needed to avoid some 
exceptional submanifolds of codimension one, which remained unknown. This ap- 
proach led to a proof of full hyperbolicity (but not yet ergodicity) for all A > 2 
and v > 2, and for generic (mi, . . . , tun, r), see [S-Sz(1999)]. Later I significantly 
simplified the arguments and made them more "dynamical", which allowed me to 
obtain full hyperbolicity for hard balls with any set of external geometric parame- 
ters (mi, . . . , mjv, r) [Sim(2002)]. Thus, the hyperbolicity has been fully established 
for all systems of hard balls on a torus. 

To upgrade the full hyperbolicity to ergodicity one needs to refine the analysis of 
the aforementioned degeneracies. For hyperbolicity, it was enough that the degen- 
eracies made a subset of codimension > 1 in the phase space. For ergodicity, one has 
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to show that its codimension is > 2. In the paper [Sim(2003)] I took the first step in 
this direction - I proved that systems of N > 2 balls on a 2D torus (i.e., v = 2) are 
ergodic for typical (generic) vectors of external parameters (mi, . . . ,rriN,r). The 
proof again involves some algebro-geo metric techniques, thus the result is restricted 
to "generic" masses. But there is a good reason to believe that systems in v > 3 
dimensions would be somewhat easier to handle, at least that was indeed the case 
in early studies. 

In this paper I will prove two theorems: 

Theorem 1. In the case v = 2, for any integer value N >2 and for every (N + 1)- 
tuple (mi, . . . , rriN,r) of the external geometric parameters the standard hard ball 

system (M^ jT ,, <^ St r \ , ^rn,r) is (fully hyperbolic and) ergodic. 



Theorem 2. In the case v > 3, for any integer value N >2 and for every (N + 1)- 
tuple (mi, . . . ,rriN,r) of the external geometric parameters the standard hard ball 

system ^M^^, |<Sl r j , ^m.rj (is fully hyperbolic and) has open ergodic compo- 
nents. 

Remark 1.1. The novelty of Theorem 1 (as compared to the results in [Sim(2003)]) 
is that it applies to each (N + 1) -tuple of external parameters (provided that the 
interior of the phase space is connected), without an exceptional zero- measured set. 

Remark 1.2. Theorem 2 speaks about exactly the same models as the result of 
[Sim(2002)], but the assertion of this new theorem is stronger than that of the 
theorem in [Sim(2002)]: Theorem 2 guarantees open ergodic components, not only 
positively measured ones. (The latter is a standard corollary of the full hyperbolicity 
in semi-dispersive billiards.) 

Remark 1.3. As it follows from the results of [C-H(1996)] and [0-W(1998)], every 
two-dimensional [y = 2), standard hard ball system (M, {S^j^R, ^i) (the models 
covered by Theorem 1) is not only ergodic, but enjoys the Bernoulli mixing property, 
as well. Likewise, in the case v > 3, every z/-dimensional, standard hard ball flow 
appears to be B-mixing on each of its open ergodic components. 



Consider the ^-dimensional [y > 2), standard, flat torus = jW as the 
vessel containing iV (> 2) hard balls (spheres) Bi, . . . , Bn with positive masses 
mi, . . . ,rriN and (just for simplicity) common radius r > 0. We always assume that 
the radius r > is not too big, so that even the interior of the arising configuration 
space Q (or, equivalently, the phase space) is connected. Denote the center of the 
ball Bi by qi G T", and let Vi = qi be the velocity of the i-th particle. We investigate 
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the uniform motion of the balls Bi, . . . , B^ inside the container T v with half a unit 
of total kinetic energy: E = - J2iLi m i\\ v i\\ 2 = ^' ^ e assume that the collisions 
between balls are perfectly elastic. Since — beside the kinetic energy E — the total 
momentum / = ^*li m i v i e ^ IS a l so a trivial first integral of the motion, we 
make the standard reduction 1 = 0. Due to the apparent translation invariance of 
the arising dynamical system, we factorize the configuration space with respect to 
uniform spatial translations as follows: (qi, ... , (/at) ~ (qi + a, . . . , qw + a) for all 
translation vectors a G The configuration space Q of the arising flow is then 

the factor torus (^T / ) N / ~) = T^" 1 ) minus the cylinders 

C id = {(qi,...,q N )e IT^- 1 ): distfe,^) < 2r} 

(1 < i < j < N) corresponding to the forbidden overlap between the z-th and j-th 
spheres. Then it is easy to see that the compound configuration point 

l<i<j<N 

moves in Q uniformly with unit speed and bounces back from the boundaries dCij 
of the cylinders Cij according to the classical law of geometric optics: the angle of 
reflection equals the angle of incidence. More precisely: the post-collision velocity 
v + can be obtained from the pre-collision velocity v~ by the orthogonal reflection 
across the tangent hyperplane of the boundary <9Q at the point of collision. Here we 
must emphasize that the phrase "orthogonal" should be understood with respect to 
the natural Riemannian metric (the kinetic energy) ||<i9|| 2 = YliLi m i\\dqi\\ 2 in the 
configuration space Q. For the normalized Liouville measure \x of the arising flow 
{S 1 } we obviously have dfi = const • dq ■ dv, where dq is the Riemannian volume 
in Q induced by the above metric, and dv is the surface measure (determined by 
the restriction of the Riemannian metric above) on the unit sphere of compound 
velocities 



{N N ~\ 

( Vl ,...,v N ) e (W) N : ^rriiVi = and ^m.ll^ll 2 = 1 I . 
i=l i=l J 

The phase space M of the flow {S f } is the unit tangent bundle Q x § d_1 of the 
configuration space Q. (We will always use the shorthand notation d = v(N — 1) 
for the dimension of the billiard table Q.) We must, however, note here that at 
the boundary <9Q of Q one has to glue together the pre-collision and post-collision 
velocities in order to form the phase space M, so M is equal to the unit tangent 
bundle Q x §> d_1 modulo this identification. 
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A bit more detailed definition of hard ball systems with arbitrary masses, as well 
as their role in the family of cylindric billiards, can be found in §4 of [S-Sz(2000)] 
and in §1 of [S-Sz(1999)]. We denote the arising flow by (M, {-S*} te R, £*)• 

In the series of articles [K-S-Sz(1989)], [K-S-Sz(1991)], [K-S-Sz(1992)j, [Sim(1992- 
I)], and [Sim(1992-II)] the authors developed a powerful, three-step strategy for 
proving the (hyperbolic) ergodicity of hard ball systems. First of all, these proofs 
are inductions on the number N of balls involved in the problem. Secondly, the 
induction step itself consists of the following three major steps: 

Step I. To prove that every non-singular (i. e. smooth) trajectory segment S^'^xq 
with a "combinatorially rich" (in a well defined sense) symbolic collision sequence is 
automatically sufficient (or, in other words, "geometrically hyperbolic", see below 
in this section), provided that the phase point xq does not belong to a countable 
union J of smooth sub-manifolds with codimension at least two. (Containing the 
exceptional phase points.) 

The exceptional set J featuring this result is negligible in our dynamical consid- 
erations — it is a so called slim set. For the basic properties of slim sets, see again 
below in this section. 

Step II. Assume the induction hypothesis, i. e. that all hard ball systems with N' 
balls (2 < N' < N) are (hyperbolic and) ergodic. Prove that there exists a slim set 
S C M with the following property: For every phase point Xq G M \ S the entire 
trajectory S R x contains at most one singularity and its symbolic collision sequence 
is combinatorially rich, just as required by the result of Step I. 

Step III. By using again the induction hypothesis, prove that almost every singular 
trajectory is sufficient in the time interval (£q,+oo), where to is the time moment 
of the singular reflection. (Here the phrase "almost every" refers to the volume 
defined by the induced Riemannian metric on the singularity manifolds.) 

We note here that the almost sure sufficiency of the singular trajectories (featur- 
ing Step III) is an essential condition for the proof of the celebrated Theorem on 
Local Ergodicity for algebraic semi-dispersive billiards proved by Balint-Chernov- 
Szasz-Toth in [B-Ch-Sz-T (2002)]. Under this assumption Theorem 4.4 of [B-Ch- 
Sz-T (2002)] states that in any algebraic semi-dispersive billiard system (i. e. in a 
system such that the smooth components of the boundary <9Q are algebraic hyper- 
surfaces) a suitable, open neighborhood Uo of any sufficient phase point xq G M 
(with at most one singularity on its trajectory) belongs to a single ergodic compo- 
nent of the billiard flow (M, {S^teR, /tx) - 

In an inductive proof of ergodicity, steps I and II together ensure that there 
exists an arc-wise connected set C C M with full measure, such that every phase 
point x G C is sufficient with at most one singularity on its trajectory. Then the 
cited Theorem on Local Ergodicity (now taking advantage of the result of Step III) 
states that for every phase point x$ G C an open neighborhood Uq of x$ belongs 
to one ergodic component of the flow. Finally, the connectedness of the set C and 
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/i(M \ C) = imply that the flow (M, {-S*} teK , n) (now with N balls) is indeed 
ergodic, and actually fully hyperbolic, as well. 

The generator subspace A it j C W /N (1 < i < j < N) of the cylinder C^j 
(describing the collisions between the z-th and j-th balls) is given by the equation 

(2.1) A itj = {(qi,...,q N )e {Wf ■■q i = q j }, 

see (4.3) in [S-Sz(2000)]. Its ortho-complement Lij C WL uN is then defined by the 
equation 

(2.2) Li j = . . . , q N ) e (W) N : q k = for k ^ and m^, + mjqj = j , 

see (4.4) in [S-Sz(2000)]. Easy calculation shows that the cylinder Cij (describing 
the overlap of the i-th and j-th balls) is indeed spherical and the radius of its base 

sphere is equal to r^j = 2r^/^^-, see §4, especially formula (4.6) in [S-Sz(2000)]. 

The structure lattice C C R" N is clearly the lattice C = (Z U ) N = Z Nu . 

Due to the presence of an additional invariant quantity / = Yl!i=i m i v ii one 
usually makes the reduction ^2f =1 rriiVi = and, correspondingly, factorizes the 
configuration space with respect to uniform spatial translations: 

(2.3) (qi,---,Qn) ~ (5i +a, . ..,gjv + a), aeT u . 

The natural, common tangent space of this reduced configuration space is 



(2.4) Z=Y_v u ...,v N ) e (W) 
supplied with the inner product 

N 

( v , v> ) = ^2mi(vi, v'i), 

i=l 

see also (4.1) and (4.2) in [S-Sz(2000)]. 

Collision graphs. Let S^'^x be a nonsingular, finite trajectory segment with 
the collisions <7i, . . . , cr n listed in time order. (Each is an unordered pair 
of different labels i,j G {1, 2, . . . , N}.) The graph Q = (V, £) with vertex set 
V = {1,2,..., iV} and set of edges £ = {ui, . . . , a n } is called the collision graph 
of the orbit segment S^ a,b ^x. For a given positive number C, the collision graph 
Q = (V,£) of the orbit segment S^ a,b ^x will be called C-rich if Q contains at least 
C connected, consecutive (i. e. following one after the other in time, according to 
the time-ordering given by the trajectory segment S^ a,b ^x) subgraphs. 



N 

E 

i=i 



rriiVi 



= 



= n a ^ = 



J<3 
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Trajectory Branches. We are going to briefly describe the discontinuity of the 
flow {<!?*} caused by a multiple collisions at time to- Assume first that the pre- 
collision velocities of the particles are given. What can we say about the possible 
post-collision velocities? Let us perturb the pre-collision phase point (at time 
to — 0) infinitesimally, so that the collisions at ~ to occur at infinitesimally different 
moments. By applying the collision laws to the arising finite sequence of collisions, 
we see that the post-collision velocities are fully determined by the time-ordered list 
of the arising collisions. Therefore, the collection of all possible time-ordered lists of 
these collisions gives rise to a finite family of continuations of the trajectory beyond 
to- They are called the trajectory branches. It is quite clear that similar statements 
can be said regarding the evolution of a trajectory through a multiple collision in 
reverse time. Furthermore, it is also obvious that for any given phase point xq G M 
there are two, w-high trees T + and 71 such that 7+ (71) describes all the possible 
continuations of the positive (negative) trajectory S^ 0,00 ^xo (S^~°° ,0 'xo). (For the 
definitions of trees and for some of their applications to billiards, cf. the beginning 
of §5 in [K-S-Sz(1992)].) It is also clear that all possible continuations (branches) 
of the whole trajectory S^ -00 ' 00 ^ can be uniquely described by all pairs (£?_, B + ) 
of infinite branches of the trees 71 and T + (£?_ C 71, B + c 7+). 

Finally, we note that the trajectory of the phase point xq has exactly two 
branches, provided that S t xo hits a singularity for a single value t = to, and the 
phase point S to xo does not lie on the intersection of more than one singularity 
manifolds. In this case we say that the trajectory of x has a "simple singularity" . 

Neutral Subspaces, Advance, and Sufficiency. Consider a nonsingular tra- 
jectory segment S^'^x. Suppose that a and b are not moments of collision. 

Definition 2.5. The neutral space Afo(S^ a,b 'x) of the trajectory segment S^ a,b 'x at 
time zero (a < < b) is defined by the following formula: 

N (S [a ' b] x) = {W eZ: 3(5 > 0) such that Va G (-5,5) 

V (S a (Q(x) + aW, V(x))) = V(S a x) and V (S b (Q(x) + aW, V(x))) = V(S b x)}. 

(Z is the common tangent space T^Q of the parallelizable manifold Q at any of its 
points g, while V(x) is the velocity component of the phase point x = (Q(x), V(x)).) 

It is known (see (3) in §3 of [S-Ch (1987)]) that N (S^x) is a linear subspace 
of Z indeed, and V(x) G Af (S [a ' b] x). The neutral space Af t (S [a ' b] x) of the segment 
S^ a ' b ^x at time t G [a, b] is defined as follows: 




It is clear that the neutral space Af t (S^ b k) can be canonically identified with 
Afo(S^ a,b ^x) by the usual identification of the tangent spaces of Q along the trajectory 
^(-00,00)^ ( gee? for instance5 § 2 of [K-S-Sz(1990)]). 
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Our next definition is that of the advance. Consider a non-singular orbit segment 
S^ a,b ^x with the symbolic collision sequence E = (o"i, . . . , cr n ), meaning that S [a ' b] x 
has exactly n collisions with <9Q, and the i-th collision (1 < i < n) takes place at the 
boundary of the cylinder C ai . For x = (Q, V) G M and W G Z, ||W|| sufficiently 
small, denote 2V(Q,V) := (<9 + W,V). 

Definition 2.6. For any 1 < k < n and t G [a, 6], £/ie advance 

a k = a(a k ): Af t (S^x)^R 

of the collision o k is the unique linear extension of the linear functional a k = 
a(a k ) defined in a sufficiently small neighborhood of the origin of Mt{S^ a,h ^x) in the 
following way: 

a(a k ){W):=t k {x)-t k {S- t T w S t x). 

Here t k = t k (x) is the time moment of the k-th collision a k on the trajectory of 
x after time t = a. The above formula and the notion of the advance functional 

a k = a(a k ) : M t (s^x) -> R 
has two important features: 

(i) If the spatial translation (Q, V) i-> (Q + W, V) (W G M t (S^ a ^x)) is carried 
out at time t, then t k changes linearly in W, and it takes place just a k (W) units of 
time earlier. (This is why it is called "advance".) 

(ii) If the considered reference time t is somewhere between t k -\ and t k , then 
the neutrality of W with respect to a k precisely means that 

W — a k (W) ■ V(x) G A ak , 

i. e. a neutral (with respect to the collision a k ) spatial translation W with the 
advance a k (W) — means that the vector W belongs to the generator space A ak 
of the cylinder C Uk . 

It is now time to bring up the basic notion of sufficiency (or, sometimes it is also 
called geometric hyperbolicity) of a trajectory (segment). This is the utmost im- 
portant necessary condition for the proof of the fundamental theorem for algebraic 
semi-dispersive billiards, i. e. Theorem 4.4 in [B-Ch-Sz-T(2002)]. 

Definition 2.7. 

(1) The nonsingular trajectory segment S^ a,b 'x (a and b are supposed not to be 
moments of collision) is said to be sufficient if and only if the dimension of 
^f t (S^ a ^x) (t G [a, b\) is minimal, i.e. dim Aft(S^x) = 1. 

(2) The trajectory segment S^'^x containing exactly one singularity (a so called 
"simple singularity", see above) is said to be sufficient if and only if both 
branches of this trajectory segment are sufficient. 
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Definition 2.8. The phase point x G M with at most one (simple) singularity is 
said to be sufficient if and only if its whole trajectory S^ -00 ' 00 ^ is sufficient, which 
means, by definition, that some of its bounded segments S^ a ' b 'x are sufficient. 

In the case of an orbit S^ 00 ' 00 ^ with a simple singularity, sufficiency means 
that both branches of S^~°° ,00 ^x are sufficient. 

No accumulation (of collisions) in finite time. By the results of Vaserstein 
[V(1979)], Galperin [G(1981)] and Burago-Ferleger-Kononenko [B-F-K(1998)], in 
any semi-dispersive billiard flow there can only be finitely many collisions in finite 
time intervals, see Theorem 1 in [B-F-K(1998)]. Thus, the dynamics is well defined 
as long as the trajectory does not hit more than one boundary components at the 
same time. 

Slim sets. We are going to summarize the basic properties of codimension-two 
subsets A of a connected, smooth manifold M with a possible boundary and corners. 
Since these subsets A are just those negligible in our dynamical discussions, we shall 
call them slim. As to a broader exposition of the issues, see [E(1978)] or §2 of [K- 
S-Sz(1991)]. 

Note that the dimension dim A of a separable metric space A is one of the three 
classical notions of topological dimension: the covering (Cech-Lebesgue), the small 
inductive (Menger-Urysohn), or the large inductive (Brouwer-Cech) dimension. As 
it is known from general topology, all of them are the same for separable metric 
spaces. 

Definition 2.9. A subset A of M is called slim if and only if A can be covered by 
a countable family of codimension-two (i. e. at least two) closed sets of ^-measure 
zero, where ji is any smooth measure on M. (Cf. Definition 2.12 of [K-S-Sz(1991)].) 

Property 2.10. The collection of all slim subsets of M is a a-ideal, that is, count- 
able unions of slim sets and arbitrary subsets of slim sets are also slim. 

Proposition 2.11. (Locality). A subset AcMis slim if and only if for every 
x E A there exists an open neighborhood U of x in M such that U D A is slim. (Cf. 
Lemma 2.14 of [K-S-Sz(1991)].) 

Property 2.12. A closed subset A C M is slim if and only if fi(A) = and 
dim A < dimM - 2. 

Property 2.13. (Integrability). If A C Mi xM 2 is a closed subset of the product 
of two smooth, connected manifolds with possible boundaries and corners, and for 
every x G Mi the set 

A x = {ye M 2 : (x,y) G A} 
is slim in M 2 , then A is slim in Mi x M 2 . 
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The following propositions characterize the codimension-one and codimension- 
two sets. 

Proposition 2.14. For any closed subset S C M the following three conditions 
are equivalent: 

(i) dim S < dimM - 2; 

(ii) intS = and for every open connected set G C M the difference set G \ S 
is also connected; 

(iii) intS = and for every point x G M and for any open neighborhood V of x 
in M there exists a smaller open neighborhood W C V of the point x such 
that for every pair of points y, z G W \ S there is a continuous curve 7 in 
the set V \ S connecting the points y and z. 

(See Theorem 1.8.13 and Problem 1.8.E of [E(1978)].) 

Proposition 2.15. For any subset S C M the condition dim S < dimM — 1 is 
equivalent to intS = 0. (See Theorem 1.8.10 of [E(1978)].) 

We recall an elementary, but important lemma (Lemma 4.15 of [K-S-Sz(1991)]). 
Let R2 be the set of phase points x G M \ <9M such that the trajectory 5'^ _00 ' 00 ' ) x 
has more than one singularities (or, its only singularity is not simple). 

Proposition 2.16. The set R2 is a countable union of codimension-two smooth 
sub-manifolds of M and, being such, is slim. 

The next lemma establishes the most important property of slim sets which 
gives us the fundamental geometric tool to connect the open ergodic components 
of billiard flows. 

Proposition 2.17. If M is connected, then the complement M\A of a slim F a set 
A C M is an arc- wise connected (Gg) set of full measure. (See Property 3 of §4.1 in 
[K-S-Sz(1989)]. The F a sets are, by definition, the countable unions of closed sets, 
while the Gg sets are the countable intersections of open sets.) 

The subsets M° and M # . Denote by M # the set of all phase points x G M for 
which the trajectory of x encounters infinitely many non-tangential collisions in both 
time directions. The trajectories of the points x G M \ M# are lines: the motion is 
linear and uniform, see the appendix of [Sz(1994)]. It is proven in lemmas A. 2.1 and 
A.2.2 of [Sz(1994)] that the closed set M \ M# is a finite union of hyperplanes. It 
is also proven in [Sz(1994)] that, locally, the two sides of a hyper-planar component 
of M \ M# can be connected by a positively measured beam of trajectories, hence, 
from the point of view of ergodicity, in this paper it is enough to show that the 
connected components of M# entirely belong to one ergodic component. This is 
what we are going to do in this paper. 

Denote by M° the set of all phase points x G M # the trajectory of which does not 
hit any singularity, and use the notation M 1 for the set of all phase points x G 
whose orbit contains exactly one, simple singularity. According to Proposition 
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2.16, the set M # \ (M° U M 1 ) is a countable union of smooth, co dimension- two 
(> 2) submanifolds of M, and, therefore, this set may be discarded in our study of 
ergodicity, please see also the properties of slim sets above. Thus, we will restrict 
our attention to the phase points x G M° U M 1 . 

The "Chernov-Sinai Ansatz" . An essential precondition for the Theorem on Lo- 
cal Ergodicity by Balint-Chernov-Szasz-Toth (Theorem 4.4 of [B-Ch-Sz-T(2002)]) 
is the so called "Chernov-Sinai Ansatz" which we are going to formulate below. 
Denote by S1Z + C <9M the set of all phase points x = (q , v ) G <9M corresponding 
to singular reflections (a tangential or a double collision at time zero) supplied with 
the post-collision (outgoing) velocity vq. It is well known that S7Z + is a compact 
cell complex with dimension 2d — 3 = dimM — 2. It is also known (see Lemma 4.1 
in [K-S-Sz(1990)], in conjunction with Proposition 2.16 above) that for ^-almost 
every phase point xq G S1Z + the forward orbit S^ 0, °°^xo does not hit any further 
singularity. (Here v is the Riemannian volume of S1Z + induced by the restriction of 
the natural Riemannian metric of M.) The Chernov-Sinai Ansatz postulates that 
for //-almost every xq G S1Z + the forward orbit S^ q, °°'Xq is sufficient (geometrically 
hyperbolic). 

The Theorem on Local Ergodicity. The Theorem on Local Ergodicity for al- 
gebraic, semi-dispersive billiards by Balint-Chernov-Szasz-Toth (Theorem 4.4 of 
[B-Ch-Sz-T(2002)]) claims the following: Let (M, {<S*} te R, ^) be a semi-dispersive 
billiard flow with the property that the smooth components of the boundary <9Q of 
the configuration space are algebraic hyper-surfaces. (The cylindric billiards auto- 
matically fulfill this algebraicity condition.) Assume - further - that the Chernov- 
Sinai Ansatz holds true, and a phase point xq G (M° U M 1 ) \ <9M is sufficient. 

Then some open neighborhood Uq C M of xq belongs to a single ergodic compo- 
nent of the flow (M, {<S*}tgR, /S). (Modulo the zero sets, of course.) 

§3. The Geometry of the Neutral Spaces 
Proof of Theorem 1 

Suppose that E = (cri, 02, . . . , u n ) is a symbolic collision sequence (for N balls) 
that is combinatorially rich in the sense of Definition 3.28 of [Sim(2002)], i. e. for 
which Corollary 3.26 of that paper applies. Assume that the exceptional, closed 
algebraic set F (of non-sufficient phase points — despite the rich symbolic sequence 
E on their orbit segment) has a non-empty, codimension-one, smooth component. 
Assume, correspondingly, that an open set Uq C M \ <9M and a number T > are 
given with the following properties: 

(1) s T (u )ndM = ®, 

(2) S^ 0,T ^x is non-singular for all x G Uq, 
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(3) E (S^x) = E = (<7i, 02, . . . , n ) for all x G t/ , 

(4) the subset 

NS (E, t/ ) = e U \ dimATo (^ [0 ' T] x) > l} 

of non- sufficient phase points is a smooth (algebraic) submanifold of Uo, 

(5) the pair of sets (Uo, NS (E, Uo)) is diffeomorphic to the standard pair 

(R 2d_1 , R 2d ~ 2 ) . 

Below we formulate (and then prove) a fundamental result, Main Lemma 3.1, 
about the scenario described in (1) — (5) above. In the sequel we will be frequently 
using the phrase "the property P(y) holds true for typical (or, sometimes called 
"generic" ) points y of a smooth, algebraic manifold N" . By this we will always 
mean that the subset 

B = {y G N\ P(y) is false} 

is contained in a proper, closed, algebraic subset of N, i. e. B can be covered by a 
finite family of proper, closed, algebraic submanifolds of N. In all such cases, the 
exceptional set B itself will be a proper, closed, algebraic subset of the considered 
ambient manifold N. 

The anticipated main result of this section is 

Main Lemma 3.1. Use all assumptions and notations from above. We claim that 
the manifold NS(E, Uo) cannot coincide (even locally) with any post-singularity 
manifold, i. e. with any open subset of S T (STZ + ) with r > 0. 

Proof. The proof will be subdivided into several lemmas. The most important 
ingredient will be the Connecting Path Formula (CPF, Lemma 2.9 of [Sim(1992)-II], 
see also Proposition 2.19 in [S-Sz(1999)]). The CPF says that for the orbit segments 
S^'^x (x G Uq, < b < T) the advance functionals oik = oc(ak) (k = l,...,m, 
E (S^ ' b h) = (01, . . . , m ) = E m , < m < n) of the collisions on x are exactly 
the solutions of a homogeneous linear system 

m 

(3.2) J2 T ^k = 0, i = l,2,...,m + P m -N, 

k=i 

where each equation in (3.2) is a //-dimensional vector equation, the coefficients 
^ik = r ifc (Si°' b ix) G W of ctfc are some well defined linear combinations of the 
relative velocities 
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(33) Afc 

Afc = vf( k ) ~ V j(k) 

(here a k = (i(k), j(k)), 1 < i{k) < j(k) < N, and v\ e W is the velocity of 
the i-th ball between the k-th and k + 1-st collisions), so that the coefficients of 
these relative velocities are given fractional linear expressions of the masses m^k) 
and rrijtk)i please see definitions 2.17-2.18, Proposition 2.19 and formula (3.27) in 
[S-Sz(1999)]. Finally, P m = P-£ m is the number of all connected components of 
the collision graph of E m = (cxi, . . . , a m ), < m < n. We note that P n = 1, 
since E = E n necessarily has a connected collision graph, thanks to the assumed 
combinatorial richness of E. 

Note that, when applying the CPF, the left-hand-side of the Connecting Path 
Formula has to be written as the relative incoming (i. e. pre-collision) velocity of 
the colliding balls multiplied by the advance of that collision. It follows immediately 
from the exposition of [Sim(1992)-II] that the equations of this sort (arising from 
all CPF's) are the only constraints on the advance functionals. The reason why 
this is true is that the fulfillment of all CPF's precisely means that the relative 
displacement (variation of position) of every pair of particles right before their new 
collision is parallel to the relative incoming velocity of these particles, and this fact 
guarantees that the variation of the newly formed relative outgoing velocity will also 
be zero, just as required by the neutrality. Our understanding is that the columns 
of (3.2) correspond to the advances «i, . . . , a m , while its rows form m + P m — N 
z/-dimensional blocks. The important feature of (3.2) is that the dimension of the 
neutral space Nq (S^x) (x G U q , < b < T, E (S^x) = E m = (<ti, . . . , a m )) is 
directly related to the rank r of (3.2): 



(3.4) dimA/'o (S [0 ' b] x^ = v(P m - 1) + m - r, 

see Proposition 3.4 in [Sim(1992)-II]. (For m = 0, i. e. when S^ 0,b ^x has no collision 
at all, the neutral space is - by definition - the <i-dimensional space of all spatial 
variations 5q of x = (q, v) at time t = 0. Of course, in that case the vector 
space Afo (S^'^x) is naturally identified with the tangent space of the parallelizable 
configuration space Q, d = dimQ = u(N — l), Po = N,r = 0,m = 0.) We note that 
the term m — r in (3.4) is the dimension of the solution set {oti, . . . , a m } of (3.2) 
which, according to Proposition 3.4 of [Sim(1992)-II], is the rank of the mapping 
$ m : No(S^x) -> M m , $ m (W) = (ai(W),...,a m (W)). As a special, though 
important, case of (3.4), we mention that if P m = 1 (i. e. the collision graph of 
x is connected), then S^ 0,b ^x is sufficient if and only if r = m—1 or, equivalently, 
if the only solution of (3.2) is the always existing trivial one corresponding to the 
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flow direction: a>i = a 2 = ■ • • = a m , see also Remark 4.5 in [S-Sz(1999)]. It follows 
from (3.4) and what has been said before that for any given positive integer D, for 
any x G Uq, and for any b (0 < b < T) the inequality 



is true if and only if the determinant detM vanishes for every square shaped sub- 
matrix M of (3.2) with the size 



vanishing of finitely many smooth algebraic functions detM. Since such a smooth 
(analytic) algebraic function detM is either identically zero on Uq or it only vanishes 
on a proper, closed, algebraic subset F of Uq, we proved 

Lemma 3.5. For every integer m (0 < m < n) there are two typical dimensions 
Di(m) and D 2 (m) such that 



for typical phase points x E NS(E, Uq). Here the threshold b m (0 < b m < T) 
denotes a selected (and then fixed) moment of time, such that S bm (U ) H <9M = 
and E (S^> b ^x) = E m = {a 1: . . . , a m ) for all x G U . 

Note that the obvious relations 






for typical phase points x G Uq, and 




(3.6) 



D 2 (m) 




hold and, consequently, we have 



(3.7) 



D 1 (m) < D 2 (m), 



m = 0, 1, . . . , n. 
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Remark 3.8. It follows immediately from the assumptions (l)-(5) (preceding Main 
Lemma 3.1) and from Lemma 3.5 that 1 = -Di(n) < D 2 {n). The inequalities 
Di(mi) > Di(rri2) are obvious for < mi < m 2 < n, % = 1, 2. 

By switching from the open set Uq to a smaller open set Uq C Uq, if necessary, 
we can assume that for the pair of sets ((Jo, NS (jjo, E^j the properties (l)-(5) 
from the beginning of this section are still true and, in addition, 



(3.9) 



D 1 {m) =dimAT {s^^x^j , 
D 2 (m) =dimAT (s [0 ' bm] yj 



for every x G Uo \ NS (Uo, E^, y G NS ^£/o, sj, < m < n. In the future the new 
pair (u Q , NS (u , E^j will also be denoted (a bit sloppily) by (U , NS (U , E)). 
Definition. Any system of homogeneous linear equations 

(3.2/a) ^r it a fc = 0(Gf), i = l,...,s 

fc=i 

will be called a permitted system if each coefficient G M." is a given linear combi- 
nation of the relative velocities and A^T (see (3.3)), such that the coefficients of 
these relative velocities are given rational expressions of degree zero of the masses 
mj(jt) and mj^) participating in the k-th collision = j(k)). A determinant 

detM will be called a permitted determinant if M is a square shaped submatrix of 
the coefficient matrix of a permitted system (3.2/a). 

We have seen that the smooth algebraic submanifold NS (Uq, E) of U can be 
defined by a single equation 



(3.10) NS (U , E) = {x G U \ detM (x) = 0} 

with a permitted determinant detMo = detMo(x). (Once a finite family of such 
determinant equations defines a smooth, co dimension-one submanifold, clearly, one 
of those equations already does the job.) 

Select a defining permitted determinant equation detM (a;) = of (3.10) with 
the property that the largest column index ko of Mq is the minimal possible one 
out of all such column indices. The index ko will be called the critical index, and 
the corresponding collision Ofc the critical collision. 

In the future we will need 
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Lemma 3.11. For any phase point x E U the following two assertions are equiv- 
alent: 

(1) There exists a point y G Uq such that 

A-(y) A-(x) A+(y) A+(x) 



\\A^y)\\ ||A^(x)||' ||A+(y)|| ||Ajf(x)| 
but at least one of the two equations 

A fc » A"» A+(y) A+ (x) 



for = 1, . . . , ko — 1, 



K»|| ||A fc -(x)||' ||A+»|| ||A+(x)| 

is false; 

(2) M, (S^o-dx) 7^ A/o (S^o^). 
Proof. One easily sees that for fixed directions 

AfOr) 



|Af(x)| 



k = 1, . . . , k - 1 



the permitted perturbations (&/i, . . . , <^v; <5f i, . . . , <5f at) of the phase point x = 
(q, v) G Uo (not changing the above unit vectors) are exactly the ones for which 
(Sq u . . . , Sq N ) G M (S^o-^x), and (Sv u . . . , Sv N ) G Af (S^o-^x) . Therefore, 
if (2) is false, then for any perturbed point y = (q + Sq; v + Sv) of x = (q, v) = 
(Sqi, . . . , SqN', Svi, . . . , Svn) (for which the equations 



A- k(v) _ A t - W A+ M _ A+(,) ^i-l,...,*-! 



IIAJWll' H A Jto)|l iiajwi 

of (1) hold true) the difference vectors Sq = (Sqi, . . . , Sqjsr), and <k> = (Sv\, . . . , 5fjv) 
belong to J\f (S [0 ^x), thus 

Af» A±» 



H A f (y)| I K»||' 

i. e. (1) is also false. 

On the other hand, if (2) is true, then one chooses an element 

Sq = (S qi ,...,Sq N ) eAf (s^ bk o-^ x ^j \Af ^ [0 ' bfc « ] x) 

(with a \ \Sq\\ small enough), and forms the perturbed phase point y = (q + Sq, v) 
for the original point x = (q,v) G Uq. It follows immediately from the relation 
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Sq G" Ao (S^^^x) that ^ (y) is not parallel to A% (x), thus finishing the proof of 
the lemma. □ 

The minimality of the critical column index k and the fact that NS (E, Uq) ^ U 
imply that 

(i) the permitted determinant detM does depend on its last column (of index 
ko), i. e. for any given collection A^ , k = 1, 2, . . . , ko — 1 of relative velocities, the 
equation detM = defines a proper (i. e. (2u — 1) -dimensional) linear subspace in 
the space of velocities (A fco , A^~ o ); 

(ii) in Lemma 3.11 the two equivalent assertions (1) and (2) are actually true for 
typical phase points x G NS (E, Uq). Otherwise the submanifold NS (E, Uq) (or, in 
other words, the validity of the equation detMo = 0) would be determined by the 
velocity history A^ , k = 1,2, . . . , ko — 1, thus violating the minimality of the index 
feo- 

We may assume, as before, that (1) and (2) of 3.11 are true not only for typical 
phase points x G NS (E, Uq), but for all points x G NS (E, Uq). 

Construction of the Transversal Curves. 

For the points x = (q, v) of NS (E, Uq) select the vectors 

Sq = Sq(x) = (5 qi ,...,5q N ) G Xq (s^o-^ \ Af (s^^x^ , 

depending analytically on x, and choose the positive numbers k = k(x), which will 
play the role of r -1 with the number r = r(x) in Main Lemma 3.1. Assign to the 
phase point x the curve 

(3.12) 7 = 7x = {n(x) = (q + tSq, v + tn5q)\ \t\ < e } 

wity a suitably small e > 0, selected uniformly for all x G NS (E, Uq). 

Remark 3.13. The reader might get confused by the fact that the linear per- 
turbation v + tn5q of the velocity leads us out of the phase space (constituted by 
phase points with a unit speed). This, however, should cause no problem in the 
proof: The entire formalism is invariant under rescaling the velocities. If somebody 
feels uncomfortably with this, he/she can rescale the velocity of ir t (x) to get a unit 
velocity vector. 

An important feature of the perturbation 7r t (x) of x is that 

A±M*)) Af(x) 
(3 ' 14) ||Af W x))||-||A±W||' k -^---' k ° I- 
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thanks to the relation Sq G Ao [S^ 0,bk o-^x) . According to the property (i) above, 
for the given values of A^(x), k = 1, 2, . . . , ko — 1, the equation detM (x) = 
defines a (2z/ — 1) -dimensional subspace 

(3.15) H x = {(A" , A+ ) G R" © R"| (A," , A) + <A+ , B) = 0} 

in the space of pairs of velocities (A^ , A^ ) eK"©K" = R 2l/ . Here the vectors 
A = A(x) G W, B = B{x) G W depend on x G Uq. Plainly, we can assume 
that this dependence is analytic, and A(x) = A(y), whenever Ao (S^°' bk o-^x) = 
Ao (S^ 0,bk o-^y). (The angled bracket denotes the standard inner product in W.) 
We observe first that, due to (3.14) above, the same vectors A, Be M" define the 
(2u — 1) -dimensional subspace H nt ^ associated with the perturbed point 7Tt(x), 
\t\ < eo- The next step in the proof of Main Lemma 3.1 is 

Lemma 3.16. For typical x G NS (E, Uo) and t G (— eo, eo) the vector 

(A- Q (n t (x)), A+M*))) 
does not belong to H x = H nt ^, i. e. 

<A-(7r t (x)), A) + (A+Mx)), B)^0. 

Proof. The function 

cf>(t) = (A- (**(*)), A) + <A+ Mx)), B) 

is clearly an analytic algebraic function of t, \t\ < eo- Therefore, the only alternative 
to an almost everywhere nonzero function (j) is that 4>{t) = 0. By using a fixed 
number (a time moment) sq slightly smaller than t ka { x ) (the time of a ko on the 
orbit of x), we express S s ° (n t (x)) in local coordinates as a function of t: 

(3.17) S SQ Mx)) = S s "(x) + ((k' 1 + s )tK5q, tnSq) , 

where Sq = Rk -iRk -2 ■ ■ ■ RiSq, and Ri is the orthogonal velocity reflection at the 
boundary <9Q at the point q (x)) . Here we used the definition (3.12) of 7Tt(x), 

and the neutrality of Sq = Sq(x) = (Sqi, . . . , Sqjy) with respect to <7i, . . . , (Jk -i- If 
we switch from (3.17) to the coordinates describing merely the relative motion of the 
i(k )-th and j(k )-th balls (i. e. we only consider q i{ko) - q j<yko) and v^ ko) - v j( fco )), 
then, in order to prove that <p{t) ^ 0, it is enough to prove 
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Sub-lemma 3.18. Let A, B be vectors in (so that at least one of them is non- 
zero), C = C2r C W be the sphere of radius 2r in W centered at the origin, Q G W 
a point lying outside of C, ^ V = Vi(k ) — Vj(k ) £*"a (velocity) vector such that 
the line / C through Q in the direction of V intersects C in two distinct points 
P = Pi and P2 (the point P being closer to Q, V pointing in the direction of q^) 5 

d = k~ 1 +sq be the distance between Q and P, and let w = K,(Sq~i(k ) ~^Qj(k )) £ ^ 
be a vector not parallel to V. For any t, \t\ < €q, let = V + tw, k C be the 
line through Q in the direction of V t , P t the point of l t D C closer to Q, and, finally, 
the reflection of Vt across the tangent hyperplane of C at Pt- We claim that 
for typical phase points x = (q, v) G NS (E, Uq) (featuring (3.17)) and for almost 
every selection of 8q — Sq(x) = (Sqi, . . . , Sqw) G A/o (5^°' 6fe o-i] x )\Ar o (S [0 ' 6fe ° ] :r) the 
smooth (analytic) algebraic function 

4>{t) = (V t , A) + (V+, B) 

is not identically zero. 

Proof. Assume that <f>(t) = 0. Let t\ < < £2 be the two values of £ G R for which 
the lines Z tl and Zt 2 are tangent to the sphere C. Since V£ = V ti (i = 1,2), according 
to our indirect assumption (V ti , A + B) = for i = 1, 2, thus the vector A + SeR" 
is orthogonal to the 2D plane IT = H(x) spanned by Vq = V = v^k ) ~ ^j(k ) an d 
w = K(Sqi(k ) ~ ^Qj(k ))- By writing A = —B + X, X ± II, the assumed identity 
= means that (V t , -B + X) + (V t + , B) = 0, i. e. 



(*) (^ + " V t , B) = 0. 

If the translated 2-plane Q + IT C W does not contain the origin, then, as it is 
easy to see, the arising difference vectors V t + — V t span the entire space W, which 
means that B = 0. So we got that either B = B(x) = for every phase point 
x = (q, v) G NS (E, Uq), or Q G II for all points x = (q, v) G NS (E, Uq) and 
every selection of 5q — 5q(x) = (Sqi, . . . , Sq^) G A/o (S^ ' h ^x) \A/" (S^' bk o] x ). In 
the first case, due to the assumtion <f>(t) = for all x G NS (E, Uq) and for every 
selection of 6q = 8q(x) = (8 qi , 5q N ) G M (S^°' b "o-^x) \ J\f (S^' b "o^ x ), we get 
that the non-zero vector A = A(x) is always orthogonal to IT, which means that the 
validity of detM (x) = (x G Uq) is already determined by the unit vectors 

k - 1 kn _ 1 

±, mi? ft- — l,...,ft-Q 1, 



l|A^(x)|| 

thus contradicting to the minimality of the critical index fc . In the second case, 
however, the equation (*) above implies that B = B(x) _L IT, so A = B + X _L IT, 
as well, and again the validity of detMo(x) = (x G Uq) would be determined by 
the unit vectors 
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k _ 1 kn _ 1 

contradicting to the minimality of k . This finishes the indirect proof of the sub- 
lemma. □ 

The successful proof of the sub-lemma completes the proof of Lemma 3.16, as 
well. □ 

The last lemma in the proof of Main Lemma 3.1 is 

Lemma 3.19. For typical phase points x G NS (E, Uq) the curve 7 X (defined in 
(3.12)) is transversal to the manifold NS (E, Uq). 

Proof. Observe, first of all, that the direction vector Sq = 5q(x) G Ao (S^°' bk o-^x) \ 
Mq (5 , [°' 6fe o]x) (defining 7 X via (3.12)) can be selected in such a way that Sq(xi) = 
5q(x2) for any pair of points xi, x^ G NS (E, Uq) sharing the same neutral spaces 
Mq (S^^x) and Mq (S[°' 6fc o-i]x). Then the angle tp{x) made by 7x and NS (E, Uq) 
at x is also an analytic function of x. The function ip(x) cannot be identically zero on 
NS (E, Uq), for in that case we would have r y x C NS (E, Uq), 4> x (t) = for each x G 
NS (E, Uq), violating Lemma 3.16. Thus ip(x) > for typical x G NS (E, Uq). □ 

Finishing the proof of Main Lemma 3.1. Let x G S T (STZ+) HNS (E, Uq). The 
proof in the paragraph of [K-S-Sz(1990)] containing the Sub-lemma 4.4 shows that 
the curve 7 X of (3.12) is tangential to the manifold S T (S1Z + ), yet it is transversal to 
NS (E, Uq) for typical x G NS (E, Uq). This proves that the co dimension-one sub- 
manifolds S T (S1Z + ) and NS (E, Uq) cannot even locally coincide, thereby finishing 
the proof of Main Lemma 3.1. □ 

Corollary: The Proof of Theorem 1. 

As it is explained in 9.3 of [Sim(2003)], it is exclusively the proof of Proposition 
3.1 (of that paper) that uses the rather involved algebraic machinery which, in 
turn, makes it necessary to drop an unspecified zero-measured set of the external 
geometric parameters (mi, . . . , mjv, r) to prove the (hyperbolic) ergodicity of the 
standard hard disk system. Our Main Lemma 3.1 above, however, proves exactly 
the statement of 3.1 of [Sim(2003)] for each hard disk system, without exceptional 
(N + l)-tuples of external geometric parameters (mi, . . . , mjv, r). Therefore, the 
results of [Sim(2003)], together with Main Lemma 3.1 of the present paper, prove 
our Theorem 1, i. e. that every standard hard disk system {y = 2) is (fully 
hyperbolic and) ergodic. □ 
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§4. What Can Be Saved in Higher Dimensions? 
Proof of Theorem 2 



Assume that v > 3, N > 2, mi,...,mjv > 0, r > 0, and consider the stan- 
dard toroidal hard ball system ^M^ )T ,, {<S^ r j, Atm,r) = (M, {S^teR, A*) with a 

connected interior of the phase space M (as always). The Theorem in [Sim(2002)] 
(more precisely, its corollary presented right afterwards) claims that every such 
system (M, {S^gR, /i) is fully hyperbolic (i. e. has nonzero relevant Lyapunov 
exponents almost everywhere), has ergodic components Ci of positive measure, and 
the restriction jS 1 * | C-i } te]R of the flow to any component Cj has the Bernoulli mixing 
property. By applying our current Main Lemma 3.1 in the case v > 3, we will prove 
the following theorem, which is actually a bit stronger than Theorem 2 (see §1) by 
also verifying the Chernov-Sinai Ansatz: 

Theorem 2/a. For every toroidal hard ball system ^M^ ir , js'i r | , Hm,r^ {N > 
2, v > 3) 

(1) the Chernov-Sinai Ansatz (see §2) holds true, and 

(2) the ergodic components Ci C M of the flow {5 , *} teM are open subsets of M 
(modulo the zero measured sets). 



We remind the reader that (2) is, in fact, a corollary of (1), the Theorem on 
Local Ergodicity for algebraic semi-dispersive billiards (Theorem 4.4 of [B-Ch-Sz- 
T(2002)]), and the full hyperbolicity of the considered standard hard ball flow (i. 
e. the theorem of [Sim(2002)]). We also point out that the essence of Theorem 
2/a is that it holds for every hard ball system (without exceptional (N + l)-tuples 
(mi, . . . , mjv, r) of the external geometric parameters), unlike the ergodicity result 
of [Sim(2004)-lj. 

Proof of Theorem 2/a by induction on N. In the case N = 2 the statements 
have been known to be true since [S-Ch(1987)], for in this case every trajectory (with 
at least one proper collision) is automatically sufficient (geometrically hyperbolic), 
and the system is actually ergodic. 

Assume now that N > 3, and Theorem 2/a has been proved for all situations with 
N' balls, 2 < N' < N. We want to prove (1) for any system ^M,^ ^, ^S 1 ^. r j- , [irn,r^ 
with N particles. 

Consider, therefore, a smooth component (an open cell) C C S7Z + of S1Z + with 
the maximum dimension 2d — 3, and let v = vc be the measure on C induced by 
the restriction of the Riemannian metric of M to C. We need to prove that for 
z/-almost every phase point x G C the forward orbit S^ ' 00 ^ is non-singular and 
hyperbolic. 
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The non-singularity of S^ 0,oo ^x (for z/-almost every x G C) is a direct consequence 
of Lemma 4.1 of [K-S-Sz(1990)]. It will turn out that, in order to prove that S^'^x 
is hyperbolic for typical points x G C, it is enough to prove 

Main Lemma 4.1. Use all of the above assumptions and notations, including 
the induction hypothesis of Theorem 2/a. We claim that for every partition n: 
{1,2,..., N} = Pi U • • • U Pfc, fc > 2, of the label set {1,2,..., N} and for every 
threshold T > the subset 



(4.2) 



i 71 = Fc.tt.t = G C| on 5^°°^ there is no proper collision 



between particles from different classes of n} 
of C has z>- measure zero. 

Proof. Since the assertion of the main lemma is a local statement, we can switch 
from C to any open subset C ^ of it, and assume that S T (C) fl <9M = 0. In the 
future the set C will be denoted by C, a bit sloppily. 

We will be studying small perturbations 5y = (5q, Sv) G T y ~M. of the points 
y = (q, v) = S T x (x G C) and, correspondingly, some orthogonal direct sum 
decompositions of the tangent spaces T y M of M at points y G <S T (C), so that 
these decompositions will "reflect" the non-interaction between particles belonging 
to different classes of n. 

By performing an easy induction on the number N — k (where k = \ir\), we see 
that it suffices to prove that 



/ 



(4.3) 



V 



7T >7T 

T'>T 



\ 



Fc,tt,t\ (J F C ,n',T' 

7T',T' 



/ 



where n' > % indicates that n' is a proper refinement of the partition n. (Of 
course, in the case k = N the union in (4.3) is empty.) According to Lemma 
4.1 of [K-S-Sz(1990)], the doubly singular phase points form a countable family of 
codimension-two submanifolds, thus typical phase points x of C do not encounter 
any further singularity on their forward orbit S( 0,oo ^x. Therefore, it is enough to 
study the points y in the subset 
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of S T (C). For the points y G F we introduce the following linear subspaces of 7^M: 
Wi = Wi{y) = {Sy = (Sq, Sv) E T y M\ Sq 3 = Svj = for j # P h 

(4.5) m 3 5( l3 = J2 m i 5V 3 = ' Yl m A V 3> = Yl m j( V 3> 5v o) = ' 

jePi jePi jeP t jePi 

(Svj) jeP . = -Bi(y) ■ (Sqj) jeP . }, 

where the non-negative operator Bi = Bi(y) is defined for the internal subsystem 
P t by (2.4) of [K-S-Sz(1990)] (the Popart 

of y plays the role of x in (2.4)) as the negative of the second fundamental form of 
the exponentially stable manifold 7exp(yi) of the Pi subsystem. The P^-subdynamics 
Sf^ ] yi of yi is a combinatorially rich forward orbit (since 

S~ T y = x^ |J F C> ',T', 

■k , T 
I 

T'>T 

yet Si yi still may be nonsufficient on a codimension-one, exceptional subman- 
ifold E. However, the tangent plane T y S T (C) of S T (C) at y contains the um- 
bilic subspace {(Sq, 6v) \ 5v = n5q} (where k = T _1 ) so, according to our Main 
Lemma 3.1 above, the manifolds S T (C) and E cannot locally coincide, i. e. 
codim (5' T (C) fl E) > 2. Thus we have proved 

Lemma 4.6. For //-almost every point y e C = S T (C) the non-negative operators 
Bi in (4.5) are actually positive definite for i = l,...,k, that is, the sub-billiard 
forward orbits Sfj^'^yi are sufficient (geometrically hyperbolic). □ 

Remark. In the case \Pi\ = 1 the linear space Wi is the null space, and the zero 
operator Bi (also acting on a null space) is considered to be positive definite, in 
accordance with the usual convention that the trivial, one-particle forward orbit 
Sj^'^Vi is considered to be sufficient. 

Next we define the "special neutral space" N n as follows: 



N n = N n (y) = {8y = (Sq, Sv) e T y M\ 8v = 0, Sq jl - Sq h = r i (v jl - v h ) 

(4.8) 

for any ji, j2 € P% with some r; 6 1, % = 1, . . . , k\. 

It is clear that the linear subspaces W\ , . . . , Wk and N n are mutually orthogonal, 
dimWi = |(|Pi| — l)v — 1| + (in the case |Pj| = 1 we have dimWi = 0), and 
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dimA^ = (fc-l)i/+#{t||P 4 |>l}, 
thus the orthogonal direct sum 

k 

(4.9) W = W(y) = N 7T ®($W i 

i=i 

has the dimension (N—l)u = d = dimQ. Clearly, the space W projects isometrically 
onto the tangent space of Q by the projection (Sq, Sv) i— > Sq, and W actually 
assumes the form 

(4.10) W = {(Sq, -B6q)\ Sq G TQ} , 

where the operator B > is the direct sum of all the positive operators Bi (i = 
1, . . . , k) and the zero operator corresponding to the space N n . Our next result is 
a fundamental transversality statement: 

Lemma 4.11. Assume that the phase point y G F (for the definition of F, please 
see (4.4)) all operators Bi = Bi(yi) are positive definite. Then the space W = W(y) 
of (4.9) is transversal to the flow-invariant hull 

J=\JS\C) 

t>o 

of C at the point y. 

Proof. Consider a normal vector ^ no = (-2? w ) £ T y TsA of the codimension- 
one submanifold J at the point y. According to (7.11) of [Sim(2003)] (or, by the 
results of §3 of [Sim(2004)-II], esp. (3.12)) we have the inequality (z, w) < 0. Since 
the operator B of (4.10) is non-negative, we have that Bw ^ z. Select a vector 
5q G TQ with the property (5q, z — Bw) ^ 0. Then the element (Sq, —BSq) of W 
is not perpendicular to no = (z, w): 

((Sq, -BSq), (z, w)) = (Sq, z) - (Sq, Bw) = (Sq, z - Bw) ^ 0. 

This finishes the proof of Lemma 4.11. □ 

So far we got that either some of the subspaces Wi = Wi(y) is transversal to J, 
or N-jt has this property. We will finish the proof of Main Lemma 4.1 by way of 
contradiction. Thus, assume that v(F + ) > 0, where 



(4.12) 



F+ = \y G F| Bi( yi ) > 0, i = l,..., A;} 
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By the transversality property explained above, either there exists an index i, for 
instance i = 1, such that W\{y) is transversal to J at y for a i/-positive set of points 
y E F+, or N n (y) is transversal to J at y for a i/-positive collection of y E F+. We 
finish the proof of 4.1 by assuming that 



The proof in the other case (when W\{y) is replaced by N- K {y) in (4.13)) is analogous. 

We have now reached the point in the proof where the celebrated "Fundamental 
Theorem for algebraic semi-dispersive billiards", i. e. Theorem 4.4 of [B-Ch-Sz- 
T(2002)] comes to play. More precisely, we will be using the so called "transversal" 
version of that theorem to construct a positively measured bundle of local stable 
manifolds, being transversal to the manifold J. Note that this "transversal" gen- 
eralization of Theorem 4.4 of [B-Ch-Sz-T(2002)] goes along the same lines as the 
transversal generalization (i. e. Theorem 3.6 in [K-S-Sz(1990)]) of Lemma 3 of [S- 
Ch(1987)]. The transversal version of Theorem 4.4 of [B-Ch-Sz-T(2002)] says that 
there exists a number e\ > such that 



Here 7} (y) = 7} (yi) denotes the local stable manifold of y\ in the sub-billiard 
dynamics {S{}. The maximum dimension of ^[^(y) (mentioned in (4.14)) is (|Pi| — 
l)u — 1 (> v — 1). Note that, necessarily, we have |Pi| > 2, otherwise the space 
Wi(y) of (4.13) could not be transversal to J. Denote the positively measured set 
featuring (4.14) simply by A. In virtue of the absolute continuity of the foliation 
7^ s ' ) ( ■ ), see [K-S(1986)], we have the inequality 



where \i = /j,^^ is the Liouville measure of the iV-ball system (M, {S^}^, /(/). 
For the phase points z G [j yeA 7i S ' ) (y) let us construct the modified forward orbit 
{Slz} t>0 by deleting all the interactions between particles belonging to different 
classes of n, thus allowing such particles to overlap. Observe that the phase points 
z of the set {J ye Ali S \y) en j°y the following property: If the labels i and j belong 
to different classes of the partition tt, then dist (q*(t), Qj(t)) > 2r — e\ for all £ > 



(4.13) 




(4.14) 




(4.15) 
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on the modified forward orbit {Slz} t>Q . However, according to Theorem 5.1 of [S- 
Sz(1999)] (or, Main Lemma 4.4 of [Sim(2002)]) this phenomenon cannot occur on 
a positively measured set of (4.15). Note that, formally, the cited results preclude 
the positively measured occurence of dist (q*(t), q*(t)) > 2r, instead of 2r — eo, but 
this little improvement of the quoted results can be achieved easily by a very small, 
basically notational, change in their proofs. The obtained contradiction completes 
the proof of Main Lemma 4.1. □ 

Finishing the proof of Theorem 2/a. We need to prove (1) of the theorem by 
using the induction hypothesis. Main Lemma 4.1 asserts that for //-almost every 
phase point x G S7Z + the following properties hold true: 

(i) S^°'°°^x is non-singular, and 

(ii) S^ ' 00 ^ contains infinitely many consecutive, connected collision graphs. 

According to Main Lemma 3.1 above (note that this is already the second use of 
that lemma in the proof of Thorem 2/a), for any x G STZ + with a combinatorially 
rich and non-singular orbit segment S^'^x, any exceptional, codimension-one man- 
ifold NS(E) (E = E (S^x) is the symbolic collision sequence of S^'^x) cannot 
even locally coincide with any smooth component of STZ + , thus the full forward 
orbit S(°'°°'x of ^-almost every x G S7Z + is non-singular and sufficient, as claimed. 
The proof of Theorem 2/a is now complete. □ 

Concluding Remark. The status of ergodicity in the case v > 3. 

What goes wrong with the proof of Theorem 1 in the case v > 3? A key step 
in proving that theorem was the use of the results of [Sim(2003)], please see 9.3 
at the end of that paper. However, the machinery of sections 4-8 of that paper 
(i. e. proving the non-existence of the so called separating, exceptional manifolds 
J, see the paragraph in §3 of [Sim(2003)] containing the properties (0)-(3) of J) 
essentially uses the assumption v = 2, as this circumstance was pointed out several 
times thruoghout those sections. Thus, the only thing that hinders the proof of 
ergodicity in the case v > 3 is a potential survival of the ominous separating, 
exceptional manifolds J. These are the only instruments capable of separating 
distinct, open ergodic components, as described in §3 of [Sim(2003)] listing the 
properties (0)-(3) of J. 
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